AD-A103  080 

UNCLASSIFIED 


NAVAL  OCEAN  SYSTEMS  CENTER  SAN  OIEGO  CA  F/G  17/7 

ON  THE  UNIQUENESS  OF  DOPPLER  TRACKING. <U) 

MAY  81  M  J  SHENSA 

NOSC/TR-684  mi 


NOSC  TR  684 


308 


Technical  Report  684 


ON  THE  UNIQUENESS  OF  QOPPLER  TRACKING 

■  ''  *** 

,  1  M  .J./Shensa 
! ,)  1  May  »81 


‘  J'j  Final  jtapart.  NovflHMB80-  Jan 


J  7"' 


r 


Prepared  for 
Naval  Sea  Systems  Command 

, •  .  /4  * 


■  -  &  $ 


AUG  i  9 


Approved  for  public  release;  distribution  unlimited 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CALIFORNIA  92152 


31  8  19  001 


NOSC/TR  684 


NAVAL  OCEAN  SYSTEMS  CENTER, SAN  DIEGO. CA  92152 


AN  ACTIVITY  OF  THE 
SL  GUILLE,  CAPT,  USN 

Commandtr 


NAVAL  MATERIAL  COMMAND 

HL  BLOOD 

Technic*!  Director 


ADMINISTRATIVE  INFORMATION 


Work  for  this  report  was  performed  from  November  1980  to  January  1981  under 
NAVSEA  funds  —  Element  62633N;  Project  F33341 ;  Task  Area  SF33341401 ;  Work  Unit 
632-WCi  7. 


Reviewed  by 
P.M.  Reeves,  Head 
Electronics  Division 


Under  authority  of 

R.H.  Hearn,  Head 

Fleet  Engineering  Department 


ACKNOWLEDGEMENT 


The  author  wishes  to  thank  P.  Reeves.  NOSC  Code  632  for  his  comments  on  the 
manuscript. 


security  Classification  of  this  page  rwhwi  Dm  em.r *<s) 


REPORT  DOCUMENTATION  PAGE 


«.  TITLE  fand  Subtitle) 

ON  THE  UNIQUENESS  OF  DOPPLER  TRACKING 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT’S  CATALOG  HUMBER 


S.  TYPE  OF  REPORT  «  PERIOD  COVERED 

Final  Report 

November  1980  -  January  1981 


«  PERFORMING  ORG.  REPORT  NUMBER 


t.  CONTRACT  OR  GRANT  NUMBERf*.) 


t.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Naval  Sea  Systems  Command 

Washington,  DC  I  is.  number  of  pages 

14 


4.  MONITORING  AGENCY  name  a  ADDRESS  (ft  dlt/ormnt  from  Controlling  0111  cm)  15.  SECURITY  CLASS,  (of  (hi*  report; 


Unclassified 


IS*.  OECLASSiriCATlON/ DOWNGRADING 

schedule 


16.  DISTRIBUTION  STATEMENT  (of  thta  Roport) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (oi  thm  abstract  entered  In  BfocJr  20,  It  dlllmrmnt  from  Report) 


19  KEY  WOROS  (Continue  on  revere*  elde  If  necasaary  mnd  Identify  by  block  numbmr) 


Doppler  tracking 
Isovelocity  receiver 


20  ABSTRACT  (Contln ue  on  rmvmram  elde  //  nmcmmmmry  mnd  identify  by  block  numbmr) 


v  Doppler  tracking  is  a  common  procedure  which  can  be  implemented  by  a  multitude  of  techniques. 

Also,  it  is  well  known  that  in  the  absence  of  a  velocity  change  in  the  receiver,  the  solution  track  is  not  unique. 
This  paper  examines  the  solution  set  of  the  Doppler  tracking  problem,  presenting  several  new  results  and 
placing  some  of  the  lore  in  a  more  rigorous  setting.  In  particular,  it  is  shown  that  (in  two  dimensions)  for 
an  isovelocity  receiver  the  solution  is  determined  up  to  a  rotation  ahd  reflection  in  the  receiver's  coordinate 
system.  If  a  single  velocity  change  is  present,  then  there  are  exactly  two  solutions.  Generalizations  to  three 
dimensions  are  also  provided.  , 


00  I  ja"M7J  1473  EDITION  of  I  NOV  «»  IS  OBSOLETE 

S/N  0102-LF -014-6601 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  of  THIS  PAGE  (Whmt  D»(»  ImwO 


SUMMARY 


Doppler  tracking  is  a  common  procedure  which  can  be  implemented  by  a  multitude 
of  techniques.  Also,  it  is  well  known  that  in  the  absence  of  a  velocity  change  in  the  receiver, 
the  solution  track  is  not  unique.  This  paper  examines  the  solution  set  of  the  Doppler 
tracking  problem,  presenting  several  new  results  and  placing  some  of  the  lore  in  a  more 
rigorous  setting.  In  particular,  it  is  shown  that  (in  two  dimensions)  for  an  isovelocity 
receiver  the  solution  is  determined  up  to  a  rotation  and  reflection  in  the  receiver’s  coordi¬ 
nate  system.  If  a  single  velocity  change  is  present,  then  there  are  exactly  two  solutions. 
Generalizations  to  three  dimensions  are  also  provided. 


CONTENTS 


1.  INTRODUCTION  .  .  .  page  1 

2.  ISOVELOCITY  RECEIVER  ...  2 

3.  MANEUVERING  RECEIVER  ...  5 

4.  EXTENSION  TO  THREE  DIMENSIONS  ...  7 

5.  CONCLUSIONS  ...  7 
APPENDIX  A  ...  9 
REFERENCE  ...  10 


ii 


1.  INTRODUCTION 


This  sliort  paper  is  concerned  with  the  following  classical  Doppler  tracking  problem: 
Given  a  source  moving  with  constant  velocity  and  emitting  a  tone  of  known  frequency  fQ. 
and  a  receiver  following  a  known  track  with  received  frequencies  fj  at  times  tj,  determine 
the  track  of  the  source  (see  Figure  la).  Assuming  the  data  are  exact,  we  wish  to  address  the 
question,  what  is  the  nature  of  the  solution  set  of  the  above  problem?  And.  specifically, 
when  is  the  solution  unique?  Although  some  of  the  results  presented  here  are  well  known 
in  the  lore  of  Doppler  tracking,  to  the  author’s  knowledge  a  rigorous  derivation  has  not 
appeared  in  the  literature.  Other  results  are  new. 

Before  proceeding,  it  is  both  mathematically  and  descriptively  convenient  to  trans¬ 
form  to  a  coordinate  system  in  which  the  receiver  is  at  rest  (Figure  lb).  For  this  to  be  an 
advantage,  we  must  assume  t hat  the  Doppler  shift  depends  only  on  the  relative  velocity 
V  =  "vs  -  v'r  of  the  source  and  the  receiver,  an  approximation  which  is  valid  as  long  as 
Iv  I  <  c  where  e  is  the  propagation  speed  of  the  signal.  In  that  case,  the  received  frequency 
is  given  by 

f:  =  f'  ( 1  +  —  cos  0  ■  )  ( 1 ) 

1  U  J 


with  8 j  the  angle  between  v  and  the  position  vector  as  illustrated  in  Figure  lb.  We  begin 
this  study  with  the  case  in  which  the  receiver's  track  is  a  single  segment  of  constant 
velocity  V  . 


Figure  1.  Example  of  source  (s)  and  receiver  (r)  geometry;  (a)  fixed  coordinate  system, 
(b)  coordinates  relative  to  receiver. 
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2.  ISOVELOCITY  RECEIVER 


A  solution  track  takes  the  form  illustrated  in  Figure  2.  It  is  immediately  clear  that 
given  one  solution,  an  arbitrary  rotation  of  the  coordinate  axes  will  yield  another;  thus,  the 
solution  set  is  infinite.  This  leads  rather  naturally  to  the  question  of  whether  the  solution  is 
unique  up  to  a  rotation  and  reflection;  i.e.,  are  the  speed  IV 1.  the  distance  to  CPA  (closest 
point  of  approach)  R,  and  tQ  the  time  of  CPA.  uniquely  determined  by  the  Doppler  data? 
Since  there  are  three  unknowns,  I VI,  R.  and  t0,  we  suspect  that  three  data  points  consisting 
of  measurement  times  tj  and  corresponding  frequencies  fj.  i  =  1,2.3.  will  determine  the 
solution(s).  (It  is  clear  that  if  the  source  is  traveling  directly  away  from  the  receiver  the 
solution  remains  indeterminate,  and  we  therefore  exclude  that  case;  i.e..  we  assume  fj  f  C 
for  i  ^  j.) 

Define  the  variables 


Figure  2.  Geometry  in  receiver’s  coordinate  system  for  isovelocity  receiver. 


Equation  (1)  becomes 


Afi  =  fi-fo  = 


‘o  lvl 


/ 


i  =  1,2.3 


,  +  qf 


and  (3)  implies 


l2-ti  _  q:~qi 

l3_ti  ^3 -  q  i 


(5) 


(6) 


Although  this  formulation  has  increased  the  number  of  unknowns  to  four  (q^  and  Ivl).  it 
shall  shortly  enable  us  to  eliminate  two. 

We  now  introduce  the  variables 


x=  — 
x  ^2 

a  qi 
y  =  tt: 


and  letting 


Af 

Tf 


">  A 

and  0  = 


Af: 

aT; 


we  have  from  (5) 


a'1  = (sgn  x) 


\A“  +  qf 


J7- 


I  VJ  +(lj 

P  -  (Sgn  y)  — 7 —  * 

VI  +qf 


<7) 


(8) 


(9) 
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or 


1  +  q|  =  or  (x“  +  qj) 

I  +  qy  =  0“  (y~  +  qf)  •  HO) 


Eliminating  qy  yields 


'>')*>  T  T  T  “>  ^ 

(1  -a-)0-y--a-(l  -0-)x“=0"-a“ 


(11) 


Also,  the  substitution  of  (7)  into  (6)  produces 


y  x  -  xy 


(12) 


or,  after  some  rearranging. 


<x- 


1  -p 


)(y  +  i 


-)  =  - 


(l  -  p)~ 


(13) 


Thus,  the  problem  has  been  reduced  to  solving  the  two  simultaneous  quadratics  (11)  and 
(13)  for  x  and  y  where  a,  0,  and  p  are  determined  from  the  data  t  j .  Afj.  i  =  1 .2,3. 

The  sign  of  q,  is  known  from  (5)  which  restricts  the  solution(s)  to  one  quadrant 
(cf.  (7)),  The  solution! s)  will  be  the  intersection(s)  of  a  section  of  a  hyperbola  (eq.  ( 13))  with 
a  hyperbola  or  ellipse  (eq.  ( 1 1)).  A  detailed  analysis  is  found  in  Appendix  A.  If  the  Afj  are 
all  of  the  same  sign,  the  solution  is  unique.  It  is  interesting  to  note  that  when  this  is  not 
true  there  exist  cases  for  which  there  are  two  (but  never  more  than  two)  solutions.  Observe 
that  if  there  are  five  data  points  available,  at  least  three  must  have  the  same  sign  and  the 
solution  of  (1 1)  and  (13)  is  unique.  Since  we  are  assuming  the  data  are  exact,  there  is  at 
least  one  solution,  and  we  conclude: 


THEOREM  I 

For  an  isovelocity  receiver,  the  source  track  is  uniquely  determined  up  to  rotations 
and  reflections  (with  respect  to  the  receiver’s  coordinate  system)  by  five  Doppler  data 
points  (satisfying  fj  f  fj  for  i  f  j).  Three  data  points  suffice  for  this  uniqueness  if  the 
Doppler  shifts  are  of  the  same  sign. 
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3.  MANEUVERING  RECEIVER 


Consider  the  ease  in  which  the  receiver’s  track  consists  of  two  segments  as  in 
Figure  1 ;  i.e..  at  some  point  in  time  there  is  a  single  change  in  velocity  from  Vf  to  vj. 
Let  AV  =  Vf  -Vp  then  the  relative  velocities  V  and  V  (7  =  Vs  -  Vr)  satisfy 


v - v 1  = -Av 


(  14) 


We  may  use  the  techniques  of  the  previous  section  to  solve  for  IVI  and  IV  I  along 
their  respective  segments  (assuming  there  are  at  least  three  data  points  for  each).  The 
triangle  formed  by  IV.  IV  I.  and  IAV  is  determined  by  its  three  sides.  Since  A"v  is  a  known 
vector,  the  orientation  of  this  triangle  is  determined  up  to  a  reflection  (see  Figure  3a).  Thus, 
given  fv  I  and  IV  '  I.  there  are  at  most,  two  possibilities  tor  Vs  =  V  +  Vf.  The  velocity  V  com¬ 
pletes  the  determination  of  the  source  track  since  it  specifies  the  orientation  in  Figure  2. 

Now.  assume  that  the  conditions  of  Theorem  I  have  been  met  for  each  of  the  two 
segments.  Then  IVI  and  IV  I  are  uniquely  determined,  and  there  are  at  most  two  solutions 
(Figure  3a)  for  the  source  track.  We  shall  demonstrate  that  there  are  at  least  two  solutions. 


Figure  3.  Illustration  of  possible  solutions  to  Av  =  v1  -  v  where  Av.  I v  I.  and  Iv'l  are  known: 
(a)  single  maneuver,  (b)  two  maneuvers. 
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THEOREM  II 


Let  the  receiver’s  track  consist  of  at  most  two  isovelocity  segments.  Then,  given  one 
solution  to  the  Doppler  tracking  problem,  there  exists  a  distinct  second  solution. 

The  proof  relies  on  the  following  lemma: 


LEMMA  1 

Given  two  vectors  v  j  and  v->.  there  exists  a  (non-unique)  vector  vQ  such  that 
v  j  -  vQ  and  v->  -  vQ  are  parallel. 

PROOF:  The  construction  is  illustrated  in  Figure  4a.  A  rigorous  proof  may  be 
supplied  by  the  reader. 

PROOF  OF  THEOREM  IT.  We  now  transform  to  a  coordinate  system  moving  with 
velocity  V0  such  that  "vr  -  vQ  and  Vr  -  v  are  parallel.  It  is  clear  from  Figure  4b  that  given 
one  solution  (solid  line),  there  exists  a  second  solution  (dotted  line)  obtained  by  a  reflec¬ 
tion  in  the  receiver’s  track. 


Figure  4.  (a)  Example  of  Lemma  1;  (b)  Illustration  of  Theorem  11. 
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Combining  these  results  with  those  of  the  previous  section,  we  have 

THEOREM  III 

Let  the  receiver's  track  consist  of  two  segments,  each  containing  at  least  five 
Doppler  points  (with  tj  f  tj  tor  1  *  j).  Then  there  exist  exactly  two  solutions  to  the  Doppler 
tracking  problem. 

Suppose  now  that  there  are  further  maneuvers.  Label  the  corresponding  velocity 
differences  with  a  subscript  i  so  that  a7(  ~~rt  -7f0  where  7  is  the  velocity  ot  the  initial 
segment.  It  is  easily  seen  from  Figure  3b  that  if  there  exist  at  least  two  differences  A7^  and 
A7t  which  are  not  parallel,  the  Doppler-tracking  solution  must  be  unique  (because  both 
~0  and  its  reflection  t dotted  line  m  Figure  3b)  in  the  triangle  containing  a7j.  cannot  match 
their  counterparts  in  the  triangle  containing  ). 

4.  EXTENSION  TO  THREE  DIMENSIONS 


The  previous  results  generalize  directly  to  the  three-dimensional  case.  Once  again 
we  transfer  to  the  receiver's  coordinates  For  ar  isovelocity  receiver,  the  source  track  and 
its  perpendicular  to  the  receiver  at  CPA  determine  a  plane.  In  that  plane  the  situation  is 
exactly  that  of  Figure  2.  only  we  now  have  spherical  symmetry.  Thus,  the  solution  is 
indeterminate  up  to  two  angles  (the  direction  of  R)  rather  than  one  angle  plus  a  reflection. 
Note  that  the  wording  in  Theorem  I  is  such  that  it  remains  valid. 

When  the  track  has  two  segments.  Figure  3a  is  still  valid  (since  7 j  and  7s  determine 
a  plane):  however,  it  may  be  rotated  about  the  vector  A  7  ( Figure  5a).  Thus,  we  find  that 
the  set  of  solutions  is  unique  up  to  a  one  parameter  set  of  rotations. 

In  the  case  of  three  segments  there  are  at  most  two  solutions,  provided  the  condi¬ 
tions  of  the  previous  section  are  satisfied:  i.e..  provided  there  exist  k  and  V  suJi  that 
AV^  and  AT^ are  not  parallel.  This  situation  is  illustrated  in  Figure  5b.  A7^  and  Av^  deter¬ 
mine  the  plane  ot  the  figure.  The  point  Q.  representing  one  solution,  does  not  necessarily  lie 
in  that  plane.  I  he  complete  solution  set  is  the  intersection  of  two  non-identical  circles,  which 
can  be  at  most  two  points.  It  7  is  coplanar  with  Av^  and  a7^  (i.e..  if  the  receiver's  track  lies 
in  a  plane)  there  is  only  one  solution.  It  the  receiver's  track  does  not  lie  in  a  plane  there  is  a 
second  solution,  namely .  the  reflection  ot  Q  in  the  plane  of  the  paper.  linallv.  for  four  or 
more  segments  the  solution  is  unique. 

5.  CONCLUSIONS 

We  have  examined  the  solution  set  of  the  Doppler  tracking  problem  for  a  fixed 
velocity  source  and  possibly  maneuvering  receiver.  The  problem  as  stated  involved  four 
unknown  parameters.  Fora  sufficient  number  of  data  points  three  of  these  parameters  are 
determined  by  one  isovelocity  receiver  track  segment.  Note  that  one  of  these  parameters 


Figure  5  (a)  Three  dimensional  case.  Tire  -.iMor  T  nuv  lie  an\  where  on  the  circle. 

(b)  fcxample  with  two  maneuvers.  The  circular  loci  do  not  lie  in  the 
same  plane  and  meet  at  the  single  point  0 

is  relative  source  speed.  A  second  segment,  ti  e.,  a  maneuver),  determines  the  fourth 
parameter  up  to  a  reflection.  In  three  dimensions,  three  segments  are  required  to  achieve 
uniqueness  up  to  a  single  reflection. 

In  addition  to  the  geometric  insight  provided,  the  analysis  may  be  used  to  find  a 
starting  point  for  the  solution  of  the  nonlinear  equations  typical  of  least-square  fits  to 
Doppler  data.  In  the  case  of  a  single  maneuver,  there  are  two  solutions,  a  piece  of  informa¬ 
tion  which  can  be  critical  in  solving  nonlinear  equations.  Finally,  we  note  that  the  construc¬ 
tion  in  Section  2  reduces  the  problem  for  exact  data  to  the  simultaneous  solution  of  two 
quadratics  which  is  much  simpler  than  the  original  problem  in  four  unknowns. 


APPENDIX  A 


Given  three  Doppler  points,  we  order  them  so  that  t^  <  ti  <  tj.  Thus. 

0  <  p  <  I.  We  shall  refer  to  the  situation  in  which  all  three  Afj  are  non-negative  as 
Case  1.  Note  that  our  conclusions  concerning  the  solution  set  in  this  case  will  remain 
valid  tor  At;  0.  i  =  1.2.3  by  symmetry.  The  remaining  possibility,  where  one  At';  differs 
in  sign  from  the  other  two.  shall  be  termed  Case  2.  Also,  from  symmetry,  we  may  assume 
without  loss  of  generality  m  Case  2  that  t->  <  Iq  <  tj. 


CASE  1 


Our  assumptions  imply  (see  (3)  and  (7|)  that 


0  <  y  <  x  <  1 


<  A- 1 ) 


and  (eq.  (8)1 


1  <  a  <  0 


( A-2) 


In  addition  (qj must  be  positive  so  that  equations  ( 10)  imply 


x  <  — 
a 


<  A-3) 


Equation  ( 1 1 )  is  therefore  a  hyperbola  with  part  of  a  single  branch  in  the  first  quadrant 
and  similarly  for  equation  ( 1 3)  (since  0  <  p  <  1 ).  In  fact  it  is  easily  shown  that  the  situa¬ 
tion  is  as  illustrated  in  Figure  A-l.  There  are  thus  at  most  two  solutions.  Note,  however, 
that  the  point  x  =  1 .  y  =  1  is  a  solution  of  both  equations.  But  this  point  does  not 
satisfy  ( A-l ).  Thus,  there  is  a  maximum  of  one  solution. 
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Figure  A- 1 .  Plot  of  hyperbolae  (11)  and  (13)  for  Case  1 . 


0  <  a  <  0  .  (A-5) 

This  gives  rise  to  three  subcases:  (a)  la  I  >  1 ,  1/31  >  1 ;  (b)  la  I  <  1 .  101  >  1 ;  and  (c)  la  I  <  1 . 

101  <  1.  Note  that  in  (b).  equation  (11)  represents  one  quadrant  of  an  ellipse.  Also,  as  in 
Case  1,  equations  ( 10)  put  bounds  on  x  and  y  since  (q  j )~  must  be  positive.  These  three 
subcases  are  pictured  in  Figure  A~2.  Note  that  there  are  at  most  two  solutions.  Examples  of 
zero,  one,  and  two  solutions  may  be  constructed  Ml. 


